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ABSTRACT 


An investigation into the general nature and properties of phase 
trajectories of third order linear feedback control systems has been 


conducted. 


The geometry of eigenvectors, eigenplanes, and the conic surface 
named the eigencone was investigated with respect to system root 
parameters. An analysis has been made of the orientation of eigenplanes 
of systems with three real roots, and of complex eigenplane geometry as 
the damping ratio C, and undamped natural frequency we were varied. 

A method is introduced for determining system roots that will locate an 
eigenplane or complex eigenplane in a pre-determined location in third 


order error space. 


Analog and digital computer programs used to obtain phase trajec- 
tories of third order linear systems are presented. The effect of root 
location in the s-plane on phase trajectories is discussed. A thorough 
analysis of the effect of initial conditions on real root and complex 
root phase trajectories is presented along with numerous photographs of 


three-dimensional phase trajectory models. 


The inter-relation of phase trajectories, eigenvectors, and eigen- 
planes is analyzed. An analytical investigation of the response of a 
complex root system to a step input was made as a first step towards 
possible application of error space geometry to discontinuous control 


systems. 
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Chapter J 
INTRODUCTION 


In recent years many papers have been written with regard to the 
analysis and design of discontinuous feedback control systems. In many 
of these papers, various schemes have been proposed for changing the 
system damping at some optimum time to produce a deadbeat response. 
However, the rigorous mathematical developments of the switching criteria 
in some of these proposals tend to obscure somewhat the overall picture. 
While this investigation does not concern itself with discontinuous 
systems, the phase trajectories of third order linear feedback control 
systems have a definite relationship to the eigenvector geometry of 
three-dimensional error space which lends itself to easy visualization 


of switching surfaces. 


The general nature of the problem was to investigate the phase 
trajectories of third order linear feedback control systems in response 
to various initial conditions. It was intended originally to devote 
the major portion of the investigation to phase trajectories of systems 
characterized by one real root and a pair of complex conjugate roots 
and only to verify Han's work (Ref. 1) regarding the trajectories of 
systems with three real roots. However, during the course of the 
investigation, it was found that the phase trajectories of both these 
types of systems have a definite relationship to the eigenvectors 
associated with the system and to certain planes and surfaces defined by 
these eigenvectors. It was found possible to describe any phase trajectory 
of any third order linear system in response to any set of initial 
conditions in terms of their relationship to the eigenvectors associated 


with the system. 


The use of eigenvectors to define surfaces and planes in phase 
Space is not a new concept but this application has been largely confined 
to mathematical definitions and manipulations The approach adopted 
during the course of this investigation was to construct various repre- 
sentative types of phase»trajectories and eigenvectors in a three-dimen- 


sional coordinate system and observe them in relation to each other. 


’ aa 





In the analysis and design of second order feedback control systems, 
in particular, those systems with non-linearities, the use of the phase 
plane has become wide-spread. It is conceivable that the use of phase 
Space could also become 4 valuable tool in the analysis and design of 


third order and higher feedback control systems. 


The basic concept of phase space is the choosing of a particular set 
of time variables of the control system as the coordinate axes of this 
Space. In most feedback control systems, the system error automatically 
seeks zero, and the phase space coordinates are normally chosen as error, 
E, error rate or velocity, E, and the (n-1) derivatives of error necessary 
for control system of order n. The behavior of the system is then 
described by the spatial coordinates E, E, E, ee Eo nese coordinates 
are the previously mentioned time variables of the control system, and 
then the system may be described mathematically by, 


og dE +Q,_, dv E = a,-d 
dt" dt" 





spent a, SE POcE = ee 


For the third order systems which this paper is concerned with, Equation 


(I-1) becomes, in more familiar form, 


E + aE + bE + cE = 0 (tT 2)) 


Thus at any time, t, there is a point in Space, called the "state" 
point, which describes the behavior of the system, and as time varies 
from t = 0 to t = ©, this state point traces a path in this three-dimen- 
Sional space, called the phase trajectory. As time, t, approaches infinity, 
the system error seeks zero, and the phase trajectory approaches the origin 


of this phase space. 


Since this three-dimensional space is based on system error and its 
derivatives, the phase space will hereafter be referred to as "error 


space". 


In previous investigations uSing phase space, for example, Bogner's 
(Ref. 2), standard matrix techniques were used to manipulate phase space 
into principal coordinate space. This coordinate system uses the eigen- 


vectors aS coordinate axes. This transformation is particularly useful 


2 





in the analysis of higher order contactor servo systems. Since this 
investigation is based primarily on physical interpretation of phase 
trajectories, this transformation to principal coordinate axes was not 


made. 


By limiting the order of the systems to be investigated to three, 
and by not transforming into principal coordinate space, physical repre- 
sentation of phase trajectories in three-dimensional error space was 
possible. To aid in this physical representation of phase trajectories, 
the coordinate planes and octants of third order error space are defined 
as in Fig. 1, and transparent lucite models of error space similar to 


Fig. 1 were constructed. 


The various phase trajectories investigated were obtained from both 
digital and analog computers. Wire models of these trajectories were 


constructed in the coordinate system represented by the plastic models. 
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Chapter 1] 
EIGENVECTORS IN THIRD ORDER ERROR SPACE 


General 

It is well known that for each real root of a third order system 
there is an eigenvector which corresponds to a phase trajectory that is 
a Straight line. (See Ref. 3.) For a system with three real and distinct 


roots, r, <r there are three eigenvectors. The eigenvector 


< 
1 ae et 


associated with the smallest root ry is referred to as the slow eigenvector 


Since the time constant of ry is larger than those of the other two roots 


and produces a slower response. Similarly, since r, has the shortest 
time constant, its associated eigenvector is termed the fast eigenvector. 
Since the time constant of ty is some intermediate value between that of 


r, and r the eigenvector associated with it is referred to as the 


1 a. 
intermediate eigenvector. Jn the case of a system defined by one real 
root and a pair of complex conjugate roots, due to the presence of only 


one real root, there will be only one eigenvector. 


Determination of Eigenvectors 


The differential equation of a linear third order feedback control 


system may be expressed in terms of the System error, E, as 
E + aE + bE + cE = 0 (Ei) 


where, in terms of the system roots i> Ty and Ty, all assumed real 


and distinct, 


a = ry tr, + r, 
b = rit, + rir, + rT 


C ey <T) 


Pee 73 


This third order differential equation may be rewritten as three 


dependent first order differential equations by first setting E = E 


7 
Then, 
E = E, =E, (TI) 
E = by RI 2. (11-2b) 
EE, *-cE, -bE, -ab, (11-2c) 





The set of equations thus obtained can be written in matrix form as 


E. = AE., i = 1,3. Thus, 
Ab al 


|e, 0 it 0 Ey 
E, =|] 0 0 he | ee EY (11-3) 
E. -Cc =-b -a E, 


The matrix of coefficients, A, has an eigenvalue A if and only if 
there exists a vector E with components ae E, and E, such that AE = E, 
where multiplication is matrix multiplication and E is considered to be 
a one column matrix. Such a vector E is called an eigenvector of the 
matrix. Eigenvectors exist satisfying AE = XE only when A satisfies the 


determinant (A - Al) = 0, where I is the identity matrix. 


=) l 0 
det (A - Al) = O -A i = 0 (11-4) 


Expanding the determinant, 


3 rae + brA +c = 0 @ikes) 


The expression obtained is the characteristic equation of the matrix, 
not to be confused with the characteristic equation of the feedback 
control system. Any root r. of this equation is an eigenvalue, and 
associated with each eigenvalue there is an eigenvector. The character- 


istic equation factors into the form 


(A, + rj), + r (A, + r,) = 0 (L1-6) 


>) 
Thus, 


To find the eigenvector associated with each root, each root is sub- 


—_ 
Stituted in turn into the matrix equation (A - NT) E. = 0. 


-). 1 0 E 

i _ 
0 -h. 1 ™ E, = 0 (iin) 
_é -b  -(ath.) pe 





-A, E. +E, = 0 (I1-8a) 


i, | 2 
“A, E, +E, = 0 (11-8b) 
~cE) “bE, - (ath. ) E, = 0 (11-8c) 


To obtain a solution, any two of the three equations may be solved in 
terms of one of the variables and an arbitrary constant K assumed for 


its value. The easiest choice here is to solve Equations (I11-8a) and 


(I1I1-8b) in terms of E, and to assume E, i 
= eee KO 
EL a ie (II -9a) 
i i 
i Rr wate) _ 
E. r. D ,K (I t=9b) 
The values obtained for i, E, and E may be checked by substitution into 


Equation I1-8c. 


The general expression for an eigenvector may then be written as 


Ey + KE, + r. KE, 


—— 
| 
ad 


e\- 


Cis 10) 
i 


where Ey E, and E, may be considered as unit vectors along their 


respective axes. In terms of E, E, E and rss where r. is any distinct 


real root of the system, this expression may be written as 


| 
* 


sora a KE - r KE (11-11) 
i | ; 


From this general expression of an eigenvector in three dimensional error 
Space, it is seen that an eigenvector has a fixed orientation which is 


dependent only upon the location of the associated real root. 


Throughout the course of the investigation when it was desired to 
determine the eigenvectors of a particular system, it was necessary to 
locate only one point on a particular eigenvector to orient it properly 
in space. From the general expression derived for an eigenvector, any 


arbitrary point on an eigenvector is given by the coordinates (E, E 12) am 


K Ke =r Kk 
amts > 5 u 
ad Se 

il 





1 om t5 = Owand 


tr, = 3.0. Here the roots are defined as positive when they are located 


in the left half of the complex s-plane. By assuming a value for K, the 


For example, assume a system with real roots at r 


coordinates of a point on each of the eigenvectors may be determined. 
Let it be assumed that K=l. The resulting point on each of the eigen- 


vectors are listed below. 


Eigenvector EB Emin 
slow (=2, Pee 
intermediate (=1 ,..15 4) 
fast (=023555 1.5) 
Note that all three eigenvectors lie in octant 2'. (Refer to Fig. 1). 
However, Suppose that K = -l. Then 
Eigenvector (E, E, E) 
slow C2 (= eee) 
intermediate Cle, =ieasy) 
fast (02333. 7-2 eo) 


These eigenvectors all lie in octant 4. A model of this example is 


shown in Fig. 2. 


If the roots of the system had been chosen as negative, i.e., 
located in the right half of the complex s-plane, then the eigenvectors 
would be located in octants 1 and 3'. Thus, for a third order system, 
all possible eigenvectors of the system are restricted to octants 


meee, 3° and 4. 


The Eigencone 


From observations of eigenvector location in three dimensional 
models, it was theorized that there is a conical surface in three dimen- 
sional error space which is the locus of all possible eigenvectors for 
a third order system. In order to prove this theory, the coordinates 
of any arbitrary point (E, E, E) on an eigenvector are written in para- 


metric form as 





E=K RE dese) (02) 





‘ Fig. 2 
igenvectorg of System with 


Three Real Roots 
he = vy, = £0 


I3= 3.0 








Eliminating the parameter K and solving for Ts 


E 


E + rE = 0 : Cr, S@ 5 Oli 13a) 
pA a ; +r, #-= (I1-13b) 
I 1 : 
1 E 
Eliminating tr. from these equations leaves 
a2 fi 
8S, 0 (11-14) 


which is the equation of a right elliptical cone. Since direct sub- 
Stitution of the parametric equations for any arbitrary point on an 
eigenvector satisfies this equation of a right elliptical cone, the 
surface of the cone is the locus of all possible eigenvectors associated 


with the roots of a third order system. 


This cone, which will be referred to as the eigencone, can be 
visualized as shown in Fig. 3 as the surface formed by an infinite 
number of eigenvectors emanating radially from the origin of the 
coordinate system. Figure 3 shows only that portion of the eigencone 
determined by the eigenvectors of roots located in the left half of 
the complex s-plane. An exploded view of the eigencone section located 
in octant 2' is shown in Fig. 4. Figure 5 shows this same section of 


the eigencone in its proper location as viewed looking into octant 4. 


From these photographs it is apparent that the apex of the eigencone 
is located at the origin of error space and is tangent to both the E and 
E axes. 1t is symmetrical about the line E = — in the E/E plane and is 
bisected by this plane. It will be noted in Fig. 5 that the surface of 
the eigencone nearest the E axis is labled SLOW while the surface nearest 
the E axis is labled FAST. The eigenvectors in these regions correspond 
to phase trajectories which have response times characteristic of these 
labels. By referring again to the parametric equations of a point on an 
eigenvector Equation (11-12, it will be noted that for r. < 1305 the 
eigenvector will be located in the SLOW region while for r. > 1.0, the 
elgenvector will be located in the FAST region. When a 1.0 the 
eigenvector is equidistant from both the E and E axes. It should be 


pointed out here that the classification of eigenvectors into slow, 


ee 


intermediate ard fast 1: not according to the labeis on the eigencone, 
but is applied only to eigenvectors associated with the roots of a given 
System. It is incorrect to arbitrarily call all eigenvectors slow if 
they are associated with r. <1, 0or tast-1f t > 1.0. For different 
Systems, the slow eigenvector of one may have a faster response than 

the fast eigenvector in the other. The labeling of the eigencone into 
slow and fast regions is intended to indicate the relative speed of 
response as related to the eigenvectors, and not to establish any sharp 


line of demarcation. 


~|]- 
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Chapter l1Tl 
EIGENPLANES IN THIRD ORDER ERROR SPACE 


General 

It has been shown (Ref. 1, 4, 5) that any two eigenvectors in error 
space define a plane, passing through the origin of error space, called 
an eigenplane (or sometimes a "hyperplane''). The eigenplane has an equation 
one order lower than the order of the characteristic equation of the 
control system. The eigenplane equation is obtained by suppressing the 
root corresponding to the pertinent eigenvector in the differential 


equation of the system. Thus, in general, 


nel wm = S 5 
Me + r-Fotl.. + AG ora Be, (III-1) 


For the third order case with real distinct roots, Yr) < Yo a Y35 


b> E + (1,473) E+ r5t. E = 0 (x, suppressed) 
Lb, = E£E + (1, +7,) E + rir, E = 0 (r, suppressed) (II1I-2) 
W, =E + (x, 41.) E + Yio E = 0 (r, Suppressed) 


are the defining equations for the three eigenplanes. 


The eigenplanes in third order error space are also defined by two 
Straight lines, the two eigenvectors corresponding to the two roots 
contained in the eigenplane equation. That is, the eigenplane defined 
by suppressing the smallest root, Ty» whose corresponding eigenvector 
is the slow eigenvector, contains, and can be defined by, the intermediate 
and fast eigenvectors. Ina like manner, the other two eigenplanes are 
defined by the combination of slow and fast eigenvectors, and by the 


combination of slow and intermediate eigenvectors 


When referring to specific eigenplanes hereafter in this paper, the 
"First" eigenplane will mean that plane containing the intermediate and 
fast eigenvectors formed by suppressing the smallest root, ry3 the 


"Second" eigenplane will mean that plane containing the slow and fast 


152 





eigenvectors formed by suppressing the root, ry 3 and the "Third" 
eigenplane will mean that plane containing the slow and intermediate 


eigenvectors formed by suppressing the largest root, r. 


Eigenplanes of Systems With Three Real Distinct Roots 


In a control system having three real distinct roots, there will, 
of course, exist three eigenplanes in error space. Since these eigen- 
planes are each defined by two eigenvectors, they should intersect the 
eigencone along lines corresponding to their two defining eigenvectors. 
This can be proved analytically by a simultaneous solution of the 


equations of the eigencone and of the eigenplane as shown below: 


= 


Equation of eigencone: E’ - EE = 0 
Equation of eigenplane: eas E + AE + BE =0 


where a ry + r5 
B = 


Sige 
Solving the equation of the eigencone for E yields: 


Z 
E = 


sb [sk 


This expression for E, when substituted into the equation of the 





eigenplane, yields: 9 
E + AE + af] =< 
E 
Or 
Z 2 


E. + AEE + BE’ = 0 


Now, assume E to have any arbitrary value, K. Hence, 


2 


E~ + AKE + BK = 0, and this expression, when solved by the quadratic 


formula, yields for E, 


- —AK + VK)? - 4BK? 
ae 


-16- 





or in terms of the roots, ry and r5 


> (nen) K + Vi?K*+ 24 Ks K-40 Ke 


Ec Z 


en Kk -gk t Vg2k?-Z kta kK 


- Z 


-wK- Ss e ck-GK + Vaik-n ky 














Z 
Thus: F = —rK, rls 
~ €? KK? K* _ -K- -K 
re ° cas ee 
oll om i = tS os r 
Hence: (E,6 = (- =, K,7 -Kr ) and ¢oo. eo) 


which are the coordinatesof any point on the two eigenvectors defining 


the eigenplane, Ww. 


If one were to visualize all three eigenplanes in error space, a 
volume would be bounded by these three eigenplanes. This volume would 
be two pyramids, in octants 4 and 2' with their apexes at the origin 
of error space and their edges corresponding to the three eigenvectors 
of the system. Fig 6 is a pictorial representation of just the 
pyramid in octant 2' defined by the three eigenvectors and bounded by 
the three eigenplanes of a system with roots, r, = 0 333, r, = 1.0, 


1 2 
and tr, = 3.0. 


The orientation of this pyramid is such that it will always lie 


within the eigencone and its size and shape will depend on the location 


-~17- 





Puramie Bounded By y Hree 


a | = 


Eigenplanes with Eigenverctors 
aT Their Intersections, 


270,333 r, © |.0 f; # 3/0 








of the system roots on the axis of reals If the system roots are all 

close together, the volume of the pyramid will be small as compared 

to the volume enclosed by the eigencone. Conversely, if ry is allowed 
= 1.0, 


to approach zero and r., is allowed to approach infinity, with r 


3 Z 
the resulting pyramid will almost completely fill the eigencone. 
If all three system roots are small, i.e., less than 1.0, the 
pyramid will be relatively small and oriented near the E axis; and, 
conversely, if all three roots are greater than 1.0, the pyramid will 


also be small, but oriented near the E axis. 


Eigenplanes of Systems Having One Real Root and a Complex Pair of Roots 


For a control system having one real] root and two conjugate complex 
roots, it can be readily seen from equations (JJI-2) which define the 
three general eigenplanes that only one eigenplane can exist for this 
system. If to and r, are of the form a + jb, the eigenplanes L, and 
b., cannot exist because of the imaginary coefficients of E and E in 
Equation (I11-1). Thus only one eigenplane, H, =E + Gee Ne =e ToT, E70) 
and, as was shown in the previous section, only one eigenvector with 


coordinates|# pat: k, +k ss] will exist in error space. This single 


i 
eigenplane existing for this type system will subsequently be 


referred tu 23s the ‘complex eigenplane". 


Closer inspection of the equation of the complex eigenplane given 


above shows that since the real root, r does not appear in the equation, 


1° 
the location of the complex eigenplane in error space is completely 


independent of the magnitude of the real root of the system. 


An investigation of the effects of varying system parameters and 
complex root location was conducted to gain more knowledge about the 
orientation of the complex eigenplane in error space In this phase of 
the investigation, the three-dimensional model was invaluable. By 
setting any one of the three coordinates of Equations (i1I1-2) equal 
to zero, the traces in the three coordinate planes of the model were 
easily computed and drawn on the transparent coordinate frame. This 
permitted easy visualization of the eigenplanes and their orientation 


with respect to each other, to the eigenvectors, and to the eigencone. 


=Ge 





The first investigation of 
the eigenplane location was 
made by keeping the damping 
ratio, C, of the system 
constant while varying the 
natural frequency of the 
System, we? keeping the 
real root of the system 
constant also. This is 
equivalent to the root 


movement shown to the right 





mo Fig. 7. 


Fig. / Root Movement ; 
wa variable, C constant 


The root locations investigated were: 


ry) = 3.0 * constant ¢ = .707 = constant 
rT) 3 oto a Ree) ee Ov (Case A) 
To 3 eer i) lO ws aie (Case B) 
T 347 7A Ora pera 0, a 22026 (Case C) 
ro 3 Seow ty oo i oe 4.242 (Case D) 


Models of the eigenplanes of these cases were built and from 
these it was seen that with C constant, varying w, causes the eigenplane 
to 'twist'’ around the eigencone. This twist can be easily seen in Fig. 8 


which shows the eigenplanes for Cases A, B, C and D respectively. 


As can also be seen from Fig. 9, views looking down the edges of the 
eigenplanes of Cases A, B, and C, the angle that the eigenplane makes with 
the eigencone appears to remain constant as ts is varied. No attempt was 
made to analytically prove this angle to be constant, nor was an attempt 
made to derive a relationship of this angle to the value of € chosen 


in the investigation. 


In the next phase of the investigation, the real root and the real 
part of the complex roots were kept constant while the imaginary part of 


the complex roots was varied. 20 





Eigen plane Movement; w, Variable, J Constant 


Case C an Or-y-t-en 
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This movement of the 
roots is shown at right in 


Pao. 10. 


The root locations in- 
vestigated were: 


f=) 3.0 * constant 





1 
ee OT 1.0 (Case B ) 
To 3 Pe Og e2 . 0 (Case F) 
13 io ey) 3.0 (Case G) 
Tp .371-0454.0 (Case H) 


Fig. 10 Root Movement; 
Im(x + jy) variable, 


Re(x * jy) constant 


In determining the equations of the traces of these complex 
eigenplanes in the coordinate planes, it was found that the trace of 
all complex eigenplanes, where they intersected the E/E plane, was 
constant. This is equivalent, in three dimensions, to a rotation of 
the complex eigenplane about a line in the E/E plane. This rotation 
of the complex eigenplane with the variation of the imaginary part of 


the complex roots is shown pictorially in Fig. 11 for Cases F, G and H. 


Without the use of the three-dimensional models, it can be seen that 
the variation of only the imaginary part of the complex roots will have 
no effect on the trace of the eigenplane in the E/E plane. The equation 
etre this trace is: 

ei Gay on) 1 ae 
and, with tr, and tr, complex conjugates, the coefficient of the E term 


in the above expression depends only on the magnitude of the real part 


of the complex roots. 


The final phase of this investigation of eigenplane geometry 
consisted of varying the system damping ratio, C€, while holding the 


natural frequency of the system, woe constant. 


Pi ee 
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Eigenplane Movement 
Im(xtjy) Variable 
Re(x 4) y) : 9 <}ak- TanTt 
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This movement of the roots is 
shown at the right in Fig. 12. 
The root locations investigated 


were: 


ry = 3.0 = constant 





w = 2.83 = constant 
n 





Fig. 12 Root Movement; 


C variable, a. constant 


T> 3 = 1.414 + j2.45 Golan (Case ]) 
> 3 22.0 +320 n= 7 07 (Case C) 
To 3 22.45 + 41.414 €¢ = .866 (Case M) 
ie = 2.735 4 jee ee 00 (Case) 


Here again, in determining the equations of the traces of these 
various complex eigenplanes, it was found that the trace of all complex 
eigenplanes, where they intersected one of the coordinate planes, was 
constant. In this instance, the constant trace was located in the 
E/E plane. This, likewise, is the same as a rotation of the complex 
eigenplane about a line in the E/E plane. This rotation of the complex 
eigenplane as the system § varies, with we held constant, is shown in 


Fig. 13, for Cases C, I,and J, respectively. 


Further analysis of the last two phases of this investigation of 
the geometry of the eigenplane has led to a generalization regarding 
the orientation of the complex eigenplane as the location of the complex 
roots of a third order system is varied. That is, as the imaginary part 
of the complex roots is varied from large values to small values, (this 
is equivalent to increasing the system €) the complex eigenplane rotates 


toward the surface of the eigencone and as € reaches 1.0, the complex 


ao 5S 
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Fig.td° 
Eigenplane Movement 
J. Variabie , Ww, Constant 


Case J. 
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eigenplane becomes tangent to the eigencone at the eigenvector corres- 
ponding to the magnitude of these two roots as they enter the axis of 
reals. At the moment of entry, the roots are real and repeated, and 
are associated with a common eigenvector. This may also be shown 
analytically by obtaining the equation of the eigenplane tangent to the 


elgencone. 


The general form of a plane tangent to a conic surface at a point is: 


(52) (E-E,) + (<3) (E-E,) + ($=) (E-E,)=O 


Bes E=E, c= C. 
Br. ei Ei, 
F-E, ea B=€. 


Ei a2 es 
where }ff(E, E, E) = E - EE (the equation of the eigencone). 


, Ee See 


Performing the partial differentiation at the point E . . 


n= EG yields 


il 
tJ 


I 
© 


AMT Oye 25) (C22) yg) 


r 


¢ oe \ 
ae C (E> Eee i) P Ay K, so » any point on an eigenvector. 
1 


qk (er S) 4 2K (E-«K) + R(E+qk) 20 
eKE + K See Ke ie +e +K" 


Combining terms and dividing through by K gives: 
iow e 26 + =e 0 
1 1G 
i 
or 


E + 2r. E E = 0 
i i 


which is the equation of an eigenplane defined by any two repeated 


real roots, r.- 


ma T= 





A special case of tangency between eigenplane and eigencone exists 
for the case of three repeated real roots. Here, the line of tangency 
corresponds to a single eigenvector which is the only eigenvector of the 


system, and is the eigenvector for all three roots. 


Extending this reasoning further helps explain the existance of the 
previously mentioned pyramid of real root eigenplanes. If a system 
gain is high enough so that the root locus consists of one real root 
and a complex pair, there will be one real root eigenvector and a complex 
eigenplane in error space. As the gain is lowered such that the complex 
roots just enter the axis of reals, the complex eigenplane is tangent 
to the eigencone at the eigenvector corresponding to the two repeated 
real roots. As the system gain is lowered further, the entering roots 
will diverge on the axis of reals and the complex eigenplane breaks 
into two real eigenplanes, now passing through the eigencone and inter- 
secting it at the real root eigenvectors. These two real root eigen- 


vectors now form the pyramid with the original real root eigenvector. 


Method of Determining Roots to Locate an Eigenplane in a Desired Location 


A method of fixing the eigenplane in error space has been developed 
whereby, given any two traces of the eigenplane on the coordinate planes, 
the required system roots can be readily determined. As has been pre- 
viously shown, only the complex conjugate roots of the system are necessary 
to define the complex eigenplane, and the real root of the system can be 
of any magnitude. The real roots, taken two at a time, will define the 


eigenplanes of a system with three real roots. 


The method will first be developed in general terms, and then a 
numerical example will be given. First, consider the differential 


equation of a third order system: 


E + aE + bE + cE = 0 (111-3) 
where: a > tr) + t5 + r, 
b = rjr, + Tor, + r4r3 


c= B53 


=99 = 





ry 
il 


cor ly 


pam al gis, See |) 


r+ 7. = 2x and r 


The equation of a typical eigenplane is given by: 


aT we 8 ae (r, + r,) E + r,r, Ewe 0 


Substitution of the above expressions for r, tr, and r,r, into w, yields: 


Z 5 Za i 


aa =E + (2x) E + ee oy) oe) 


This expression will give the following traces of al in the three 


coordinate planes: 


rn E/E plane) 2xE age + voy Ee = (III-4a) 

(In E/E plane) E + 2xE an ©) (II1-4b) 
oe ee 2 

(In E/E plane) E Pa e ty ) E=0 (III-4c) 


Rewriting the above trace equations in terms of the angle at the origin 


yields: 
E x? + y 
a =-lanx = De (III-5a) 
EF , 
>) hh ON) = s 
Hs fe) me (11I-5b) 
Mons ; 
es =-tan Y = y? (I11-5c) 


=). 0= 





Where a, 8, and vy are defined in Fig. 14 





Fig. 14. Angles a, B, and y on the coordinate planes. 


By determining any two of the three angles, a, B, and y, and 
substituting these angles into the corresponding tangent equations(III-5), 
a system of two equations in the two unknowns, x and y,are established. 

On solving for x and y, t, and r, are determined. If y is determined to 
be complex, substitution into the expression ro 3 = x Holy )yselds To 

and r, as real roots. If y is real, then T, and r, are complex conjugate 
roots. 


Numerical Example for Complex Eigenplane 


Assume that a complex eigenplane is required to be located such 


that its trace in the E/E plane is defined by @ = Sea and in the 


E/E plane by 8B = 60°. 


Substituting these into Equations (II1i-5a) and (III-5b) yields: 





Vee 
tan 45 = Zon ei 
Fook (I1I-6a) 
° il _ \ i. 
tan GO = eae 2 
(III-6b) 


Solution of (II1-6b) directly determines x; 
ee aati 
ray 
2306 





and this, when substituted into (I]I-6a) yields: 


(4.288) + y* = 42 (4.288) 


— awe 6 =.0005m=—2. 4927 


or y= Y.4927 = .701 


Therefore the required roots ty and r, are: 


ry 3 =2250 a ae Oe 


Numerical Example for Real Root Eigenplane 


Assume that an eigenplane is required to be located such that its 
trace in the E/E plane is defined by a = 15°, and in the E/E plane by 
= 45°. 


Substituting these values into Equations (III-5a) and (III-5c) yields: 





2 2 
so os = 268 
Fatt aX 


tan 45+ ey = 


Substitution into the upper equation yields x directly: 


Ee 268 
LIS 
> JS, See) 
2 (-268) 536 


Gecs) 4" = | 


yo? }— 3.475 
4 = +\ -2.4¢75 
15S 


y = a | 


ey ee 


ae) 


si 


Pecs ti ( jh 573) 
pee 
r, = 1.865 -j(j573) 


- 3.439 


Note that since ty and + 


3 are complex conjugate roots, only one value 
of y, either positive or negative, need be substituted into r 


23 ee 


290% 





Chapter IV 
PHASE TRAJECTORIES OF THIRD ORDER LINEAR SYSTEMS 


General 

It is well known that the time solution to the differential 
equation of a linear third order feedback control system can be repre- 
sented by a continuous curve in three-dimensional error space. If the 


differential equation of a third order system, 
E + aE + bE + cE = 0 (IV-1) 


is solved by Laplace transform methods, the time solution obtained is: 


E (t) _ Eee + (y,+1,)E, + Ge Vek. -v t E. + (G+ E+ (y, ae =a t 





e =f 
a tema Creaaee 
E um). dete E., ~~, 0 
ae feo Sane) é 
See le) (Iv-2) 


where E> ae and E are the initial conditions of error, error velocity, 
and error acceleration respectively. Successive differentiation of 


Equation (IV-2) with respect to time yields similar expressions for E(u) 


and E(t). 


Substituting the values of initial conditions ES? EO and oe system 
root values, and time into Equation (1V-2) and its derivatives yields the 
coordinates of points which describe the "state", i.e., the error, error 
velocity, and error acceleration of the system at that particular time. 


The locus of all these points is the phase trajectory of the system. 


The calculation of trajectories by direct substitution into 
Equation (IV-2) and its derivatives would be a laborious process. This 
simple but time consuming method of solving for phase trajectories 


suggests solution by means of electronic computers. 


aan 





In the initial stages of this investigation, solution of the 
differential equation (IV-1) was accomplished utilizing the Control Data 
Corporation 1604 high speed solid state general purpose digital computer. 
The digital computer program was written in FORTRAN language and is 
given in Appendix A along with its flow diagram. This FORTRAN program 
is based on the Runge-Kutta method of numerical integration. In order 
to obtain the desired accuracy of solution, state points were computed 
at intervals of 0.01 sec. but were printed out at intervals of 0.05 sec. 
A sample print out of the digital solution is also given in Appendix A. 
The program was allowed to run for a total problem time of 10.0 sec., 
and in practically all cases, at the end of this time, the values of 
Ss E, and E were small enough to consider that the trajectory had reached 
the origin of error space. In the beginning of the investigation, it had 
been anticipated that a digital to x-y plotter would be used to obtain 
the trajectories in graphical form. However, this plotting equipment did 
not become available, and it was necessary to hand plot all phase 


trajectories. 


In the latter stages of the investigation, in order to obtain 
trajectories more rapidly, a Donner 3100 analog computer was used in 
conjunction with an x-y servo plotter. The flow chart and scaling data 


for the analog computer solution are given in Appendix B. 


In both the digital and analog solutions, plots of the trajectories 
in the three coordinate planes, E/E, E/E, and E/E were made and then 
transferred to the corresponding planes of the plexiglass models. Using 
these projections, the three-dimensional models of the phase trajectories 
were constructed. Samples of the three projections of a typical phase 
trajectory of a three real root system, as obtained from the analog 


computer, are given in Figs. 15 through 17. 


Brrect Of Root Location on Phase Trajectories 


For a system defined by a given set of roots, the phase trajectories 
in response to various initial conditions are constrained by the eigen- 
vectors and eigenplanes associated with the defining roots. Since it has 


been shown previously that eigenvector and eigenplane orientation are 
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dependent only on root location, any change in root location will cause 
a reorientation of eigenvectors and eigenplanes and so define a different 
System. The phase trajectories of this system will be altered from 


that of the original system in a manner dictated by the new orientation. 


For a system defined by three real roots, it has been previously 
shown that any change in root location will effect the orientation and 
size of the pyramid bounded by the three eigenplanes. The phase tra jec- 
tories will exhibit the same general characteristics in response to initial 
conditions as discussed in the following sections but will be dependent 


on the particular orientation and size of the pyramid. 


AS was previously indicated, a system defined by one real root and 
two complex conjugate roots has only one eigenvector and one eigenplane. 
Changing the location of the real root changes the time response of the 
system by moving the eigenvector along the surface of the eigencone 
towards either a slower or faster trajectory. When the complex roots 
are changed, the eigenplane is reoriented and the phase trajectories are 
changed with regard to their E, E, and E overshoots. In addition to 
eigenplane reorientation, if the complex roots are changed So as to 
produce a change in damping ratio, C, the general oscillatory nature of 


the trajectories is increased or decreased. 


Berect of Initial Conditions on Phase Trajectories 


One of the primary purposes of this investigation was to study the 
effect of varying initial conditions on the phase trajectories of third 
order linear systems. In particular, three general areas for the 
location of these initial conditions were selected: 

1) Initial conditions on an eigenvector. 

2) Initial conditions in an eigenplane. 

3) Initial conditions not in an eigenplane. 


Each of these will be discussed in the following paragraphs. 


l. Initial Conditions on an Eigenvector. 
If initial conditions corresponding to any point located on an eigen- 
vector are substituted into Equation (1V-2), two of the three terms on 


the right side of the equation vanish. AS an example, let: 


492 





», E.= -K, and E = Kr 
O 


Substitution into Equation (IV-2) yields: 


Y, V3 
ft er am -yt 


From which 


o¢ Mie 
F (t) = ke and F (4)= K i 


Hence it can be seen that at any time, t, the state point of the 


: K 
system has the same coordinates, pa -K, Ke, | » aS any arbitrary point 
1 

on an eigenvector and the trajectory is a straight line coinciding with 


the eigenvector. Figures 18 through 20 are typical analog solutions of 


a trajectory with initial conditions on a system eigenvector. 


This statement that a phase trajectory having initial conditions on 
an eigenvector will remain on the eigenvector holds for both a system 
with three real roots, and a system with one real root and two complex 


conjugate roots. This is apparent from Equation (IV-2) which holds for 


both systems. 


As an interesting sidelight to this property of phase trajectories 
with initial conditions on eigenvectors, an unstable system with complex 
roots in the right half plane and its real root in the left half plane 
was simulated on the analog computer. This system was given a set of 


initial conditions on an eigenvector, and as the solution was generated, 
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the phase trajectory coincided with the eigenvector until just before 
it reached the origin. Here, the instability of the system caused the 
development of diverging oscillations about the origin. Figure 21 is 


the projection of this trajectory on the E/E plane. 


2. Initial Conditions in an Eigenplane. 

If initial conditions corresponding to any point in an eigenplane 
or a complex eigenplane are substituted into Equation (IV-2), Han has 
shown (Ref. 1, pp. 26-36) that the resulting phase trajectory will 


always remain in the plane. 


A phase trajectory starting in a complex eigenplane will be an 
oscillatory curve that spirals into the origin of error space, always 
remaining in the complex eigenplane. Figure 22 shows the general 
oscillatory nature of the trajectories of four different systems having 


one real root and two complex conjugate roots. The systens illustrated 


are: 
Case L vr, 20.5 ro 3 7124 33 (CG #®. 3) 
Case B r, = 3.0 rr 3* Tee Corr 7) 
Case C r, = 3.0 To 37 cot Ce 3,70) 
Case Ar, * 3.0 ty 3 My, See 8 CG 707) 


Figure 23 is a view of the same model, but looking down the edge of the 
complex eigenplane of Case A. Figure 24 shows the phase trajectories of 
Case L only, looking into octant 1. Figure 25 is a view of the same 
model of Case L but looking down the edge of the complex eigenplane. The 
dashed line in Figs. 24 and 25 represents the eigenvector corresponding 


to the real root, ry = 0.5, of Case L. 


In a system having three real and distinct roots the phase trajec- 
tories will also remain in an eigenplane if the initial conditions are 
located in this eigenplane. However, the trajectories will not have the 
oscillatory nature as in Figs. 22 through 25, but will be curves that 
become asymptotic to the slower of the two eigenvectors which define 


the eigenplane. 
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Figure 26 shows a model constructed of a rumber of phase trajec- 


tories of a system with three real roots. 1, = 0.5, r, =]1 0, 


] 2 
4 = 2.0, all having initial conditions in the Third eigenplane, defined 
by the slow and intermediate eigenvectors. Note here that two of the 
phase trajectories, (A) and (B), are straight lines in octants 4 and 2', 
indicating that the initial conditions for these trajectories were 
chosen on the slow eigenvector. This tigure readily illustrates the 
above statement that the phase trajectories become asymptotic to the 


slower of the two defining eigenvectors. 


Figure 27 is another view of the same model but looking down the 
edge of the Third eigenplane, and showing that the phase trajectories 


all remain in the eigenplane. 


Figure 28 shows phase trajectories of a system having three real 


Bf. 333, 23 Ot oO ein this figure, Trajectory.) 


1 Zz S| 
has initial conditions in the First eigenplane and becomes asymptotic 


roots, r 


to the intermediate eigenvector, Trajectory (B) has initial conditions 
in the Second eigenplane and becomes asymptotic to the slow eigenvector, 
and Trajectory (C) has initial conditions in the Third eigenplane and 


also becomes aSymptotic to the slow eigenvector. 


3. Initial Conditions not in an Eigenplane. 

Han and Thaler have stated (Ref. 3) that eigenplanes "subdivide 
the (error) space into regions so that no phase trajectories can pass out 
of one region and into another. Thus the (eigenplanes) act as boundaries 


which funnel the phase trajectories into the origin". 


For a system having three real and distinct roots, this statement 
can be modified to read that "the pyramid of eigenplanes subdivides 
error space into regions so that phase trajectories starting inside 
the pyramid will remain inside and those phase trajectories originating 
outside the pyramid will never enter the pyramid of eigenplanes. Thus 
the pyramid acts as a boundary which funnels the phase trajectories into 


the origin". 


It is important to note that not only does the pyramid funnel the 


phase trajectories into the origin, but it funnels them in asymptotic 
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to the slow eigenvector, no matter how close the initial conditions 
might be to the fast or intermediate eigenvectors. Figure 29 shows 

two phase trajectories of a three real root system with roots: 

ty ae), 333, t, Pele (.. tr, = 3.0, whose initial conditions are inside 

the pyramid of eigenplanes. It can readily be seen that both trajec- 
tories approach the slow eigenvector asymptotically even though (A) has 
its initial conditions near the intermediate eigenvector and (B) has 


its initial conditions approximately halfway between the intermediate and 


fast eigenvectors. 


Figure 30 is an illustration of a number of phase trajectories 
whose initial conditions are all outside the pyramid of eigenplanes. 
Here again it can be seen that all phase trajectories originating 
outside the eigenplanes approach the origin of error space asymptotic to 
the slow eigenvector. Three of these trajectories exhibit interesting 
properties. Trajectory (A) originates in octant 4 near the fast eigen- 
vector, but does not approach the origin along the slow eigenvector in 
the same octant. Instead it penetrates into octants 3 and 2, then 
finally approaches the origin asymptotic to the slow eigenvector in 
octant 2'. Trajectory (B) originates in octant 4', but instead of going 
to the slow eigenvector in the adjacent octant 4, it journeys through 
octants 4, 3, 2, and finally approaches the origin along the slow eigen- 
vector in octant 2'. Trajectory (C) originates near the intermediate 
eigenvector in octant 4, approaches this intermediate eigenvector, 
parallels it for a short time, and then finally becomes asymptotic to 


the slow eigenvector very near the origin in octant 4. 


These unusual excursions through error space of trajectories whose 
initial conditions are outside the pyramid of eigenplanes are results 
of the basic property of eigenplanes, i.e., eigenplanes act as boundaries 
on phase trajectories. When the initial conditions are between any two 
eigenplanes, the phase trajectories must always remain between them until 
the trajectories reach the slow eigenvector which exists between or on 


these bounding eigenplanes. 
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In the investigation of the effect of initial conditions not in an 
eigenplane, more attention was focused on those systems having one real 
root and two complex conjugate roots. A number of interesting results 


were obtained. 


The first complex root system which was studied had roots of: 
1 sao. 0 , ro 3 


having a value for € of 0.05. As might be expected, the phase tra jec- 


a = .25 * j>.0." This system was very lightly damped, 

tory was highly oscillatory. This can readily be seen in Fig. 31. 

In this particular three-dimensional model, the scaling of E, E, and E 
values of points on the phase trajectory were modified in order to 
obtain a model capable of analysis. Here, for one coordinate axis unit, 
Eo, 2, E = 4, and E = 1. When a random initial CondLei One roint (Aye 
not in the complex eigenplane, was applied to this system, the phase 
trajectory very quickly became asymptotic to the complex eigenplane. 
Figure 32, a view looking down the E axis, shows that the phase trajec- 
tory has become asymptotic to the complex eigenplane after less than one 
oscillation. Here it should be noted that this model pictured in Figs. 31 
and 32 does not show the last five oscillations of the complete tra jec- 
tory, and they would continue from Point (B) and eventually reach the 


origin of error Space. 


The next complex root system studied had roots at ry =O. 
rT) 3 = 1.0 + j3.0. The system was lightly damped, having ¢€ = 0.3. 
The general nature of this system's oscillatory trajectories was seen 
previously in Fig. 24 In this investigation, a variety of initial 


conditions not in the complex eigenplane were applied and analyzed. 


The first initial condition applied was a step input In this 
instance, the step had a magnitude of 4 0 radians. Curve (A) of 
Fig. 33 is the resulting phase trajectory. Also shown in Fig. 33 is 
the real root eigenvector represented by the dashed line, and two 
trajectories having initial conditions in the complex eigenplane, 
Curves (B) and (B'). Figure 33 indicates that for this system, the 
phase trajectory starting outside the complex eigenplane does not 


approach the origin of error space asymptotic to the complex eigenplane, 
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but instead appears to become asymptotic to the system eigenvector 

Figure 34 is another view of the same model, but looking up the eigen- 
vector in octant 4. Note here the similarity of the general shape of the 
phase trajectories having initial conditions in and out of the complex 


eigenplane. 


The next initial conditions applied to this system was an initial 
velocity of magnitude 2.0 rad/sec . Curve (A) of Fig. 35 is the resulting 
phase trajectory. Also shown are the system eigenvector and the two 
trajectories in the complex eigenplane as before. Fig. 36 is a view 
of the same model looking down the edge of the complex eigenplane. 
Although it is difficult to tell from Fig. 36, the phase trajectory did 
not immediately become asymptotic to the complex eigenplane, as it did 
in the preceding system having € = 0.5. Figure 37 is another view of 
the same model looking down the eigenvector in octant 2'. Note again the 
similar curvature of the phase trajectories having initial conditions 


in and out of the complex eigenplane. 


Next, an initial condition very near the eigenvector was applied to 
the system. As can be seen from Fig. 38, the phase trajectory quickly 


became asymptotic to the eigenvector instead of to the complex eigenplane. 


Finally an initial condition that was equidistant from the eigen- 
vector and the complex eigenplane was applied to the system. Figure 39 
shows the resulting phase trajectory as Curve (A) and the two trajectories 
in the eigenplane, curves (B) and (B') as well as the system eigenvector. 
As can be seen in Fig. 39, the phase trajectory again becomes asymptotic, 
not to the complex eigenplane, but to the eigenvector, developing the 
characteristic spiral around the eigenvector as was seen previously in 
Figs. 33 and 38. Figure 40 is a view of this model looking up the 
eigenvector in octant 4. The similarity of the phase trajectories 
with initial conditions in and out of the complex eigenplane noted here 
and in previously discussed views looking up the eigenvector indicated 
that there must be some relation between the phase trajectories and the 


eigenvector. This led to the next phase of the investigation. 
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However, before proceeding further, the results of this phase of 
the investigation indicate that for systems having one real root and two 
complex roots, the phase trajectories of the systems that are very lightly 
damped tend to become asymptotic to the complex eigenplane. For systems 
with moderate to heavy damping, the trajectories become asymptotic to 


the system eigenvector as they approach the origin of error space. 


Interrelation of Phase Trajectories, Eigenplanes, and Eigenvectors 


In the study of phase trajectories of systems having one real and 
two complex conjugate roots, it was noted that a striking similarity 
existed between trajectories with initial conditions inside and outside 
the complex eigenplane when the trajectories were viewed looking along 
the system eigenvector, as in Fig. 40. Further investigation led to 
the following development of an interrelation among phase trajectories, 


eigenplanes, and eigenvectors. 


It can be seen from the mathematics of Chapters I] and III, and 
also in Figs. 24 and 25, that the system eigenvector pierces the complex 
eigenplane at the origin of error space. Mathematically, the complex 
eigenplane can be translated along the eigenvector, keeping the origin 
of the complex eigenplane coincident with the eigenvector, and keeping 
the complex eigenplane oriented parallel to its original position. A 
set of initial conditions can be applied to the system corresponding to 
the point on the eigenvector containing the origin of the translated 
eigenplane. If the complex eigenplane is then allowed to move along 
the eigenvector at the same velocity as the system state point, this 
state point will always coincide with the origin of the complex eigen- 
plane as these points move toward the origin of error space. Extending 
this further, let the complex eigenplane be translated again along the 
eigenvector, and impose a set of initial conditions on the system such 
that the initial state point lies in the translated eigenplane some 
distance from the translated origin. I1f the complex eigenplane is 
allowed to move toward the origin along the eigenvector as before, the 
state point now traces the phase trajectory in three-dimensional error 


Space. What has really happened is that the state point has never left 


she 





the translated complex eigenplane The trace 1t has made in the complex 
eigenplane is exactly the same as the phase trajectory described by the 
state point had the initial conditions been applied in the original 
untranslated complex eigenplane at the same distance from the origin as 
the present initial conditions were from the eigenvector in the trans- 
lated plane. Hence it can be seen that any phase trajectory is the 
result of the superposition of two separate sets of initial conditions, 


one set on the eigenvector and the other set in the complex eigenplane. 


In a system with three real roots, the same results are obtained if 
an eigenplane is translated along the eigenvector not defining the eigen- 


plane. 


Hence, any phase trajectory in three-dimensional error space can 
be described as the superposition of the two trajectories resulting from 
initial conditions on an eigenvector, and initial conditions in either a 


complex or real root eigenplane. 


This hypothesis was then checked using the analog computer set-up 
as given in Appendix C. In Fig. C-1, the upper system was given initial 
conditions corresponding to a point lying on the eigenvector. The lower, 
or "primed" system was given a set of initial conditions corresponding to 
a point offset from the eigenvector by a distance equivalent to the initial 
conditions in the complex eigenplane. The difference of the two systems 
was then plotted by x-y recorder and compared with the phase trajectory 


in the complex eigenplane. 


Figure 41 is the E/E projection of the phase trajectory in the 
complex eigenplane and Fig. 42 is the E/E projection of the phase tra- 
jectory of the "difference" system. The validity of the above hypothesis 
can be seen by a comparison of the identical projections of the phase 


trajectories of Figs. 41 and 42. 


Complex Root Phase Trajectories in Response to a Step Input 


In Fig. 33, the phase trajectory of a system with one real root and 
a pair of complex conjugate roots is shown for a step input. By noting 
the tight spiral of the trajectory around the eigenvector, it must be 


concluded that for this particular case, the trajectory intersects the 
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eigencone. Since trajectories of this type could prove of great impor- 
tance in switching applications, it was decided to investigate additional 
trajectories with step inputs to see if a definite pattern could be 
established as to if, where, and when the trajectories intersect the 


Sigencone . 


The three-dimensional models which had proved invaluable heretofore 
could not be constructed precisely enough tc provide anything more than 
a general description of the trajectory However, from the models it was 
established that all complex trajectories do not intersect the eigencone. 
A trajectory of this type is shown in Fig. 31. Also, a study of Figs. 5 
and 33, indicates that if a trajectory does intersect the eigencone, it 


will occur at a faster eigenvector than that of the system 


After having exhausted the possibilities of the three-dimensional 
models, an analytical approach was adopted. If the initial conditions 
of a phase trajectory are chosen such that E. = Ey =OOne Chey Lait len) 
condition EY is an initial step and can be considered as a step input. 
Introducing these initial conditions into Equation (IV-2): 


ees a Pg 6 


Sena (co eC eee ee rm 


(on \(tg- 1) (Ni-ta)(t3-%) — 9 (Iv-3) 


Taking the first and second derivatives of the above expression: 


= “A A By cia 
F(t): E 2 Giaty 6) Cee eee BC tae nee ce : 
7 “| (-n)(G-h) (ri-7)(G-2) (11 "3)( 12-13) (Iv-4) 
= cS -rt - -rt ae ant 
Oa na ee Peete (ie ca 
(ana) (W-n)(5-8) (r,-¥3)(ra- 3) (IV-5) 
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Substituting the above expressions for E(t), E(t), and E(t) into 


Equation (11-14), the equation of the eigencone: 


, “(ne aje aie yt -(r+m)t 


aah") (t2-1)( te) (F.-12)("-3) 





(IV-6) 


This equation satisfies the conditions that exist when and if the phase 
trajectory in response to a step input intersects the eigencone. It 

can be seen from this equation that if an intersection occurs, the time 
of intersection is independent of the size of the step input. Further- 
more, on obtaining a time of intersection, the equation becomes linear 
and indicates that all points of intersection are along a straight line 


on the surface of the eigencone. This line, by defimition, must be 


along the same eigenvector. 


In order to obtain a solution to Equation (IV-6), the equation is 


simplified as shown in Appendix D to obtain the transcendental equation: 


Cos Bt +o Sin Bt =e (a-a) t (IV-7) 
where: a = real root 
@ = real part of complex root ‘ 


B = imaginary part of complex root 
De 2 "af 
ga_-o = Bo 3 


Ba 


If a time of intersection exists other than the trivial case 
of t = Q, it can be obtained from a graphical solution of the above 
equation. If a solution exists, in all except special cases, it is 
possible to obtain more than one solution if the plot is continued 
beyond the first intersection. Figure 33 is a trajectory which would 
produce multiple intersections as noted by the tight spiral around 


the eigenvector. 





In order to avoid plotting Equation (IV-7) for systems with no 
intersections, R. L. Ashford (Ref. 6) has devised a criteria whereby 
it is possible to determine whether or not a solution exists. In 
essence, it amounts to expanding both sides of Equation (IV-7) ina 
Maclaurin series and taking the coefficient of the first power term in 
order to determine the slope of each curve. The slope of the curve for 
the left side of the equation is given by oO8 while that of the right is 
(a-@). Ashford states that when (a-@%) is less than zero, there will 
always be an infinite number of solutions, and when (a-@) is greater 
than zero, a non-trivial solution will exist provided OB is greater 
than (a-q@). Ashford's criteria, then, states that a phase trajectory 
in response to a step input will intersect the eigencone whenever: 
1. (a-a~) <0 
2.  (a-@) > 0 and of > (a-qa) 
If, through Ashford's criteria, an intersection is indicated, a 
graphical solution of Equation (IV-7) will determine the time of inter- 
section. By substituting this value of time into Equations (IV-3), 


(IV-4), and (I1V-5), the point of intersection of the trajectory with 


the eigencone is determined and hence, the eigenvector. 
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Chapter V 
CONCLUSIONS 


From the results of this investigation, the following conclusions 
can be made regarding the eigenvector geometry and phase trajectories of 
third order linear feedback control systems: 

1) The orientation of system eigenvectors in third order error 
space if a function of the location of the roots of the system on the 
complex s-plane. 

2) The locus of all eigenvectors associated with a third order 
system is a conic Surface describing two right elliptical cones located 
in octants 1, 2, 3' and 4' with apexes at the origin of error Space. 

3) An eigenplane can be fixed at a desired orientation in error 
Space by proper choice of system roots. 

4) The orientation of the complex eigenplane in error space is 
completely independent of the real root of the system and depends only 
on the location of the complex conjugate roots on the complex opens 

5) The complex eigenplane never intersects the eigencone except at 
the origin, however, it can become tangent, and its point of tangency is 
along the eigenvector associated with a pair of repeated real roots. 

6) The three eigenplanes of a system having three real and distinct 
roots form a pyramid interior to the eigencone. This pyramid serves to 
funnel phase trajectories into the origin of error space. 

7) Phase trajectories of systems with three real roots having 
initial conditions in an eigenplane will remain in this plane and will 
approach the origin of error space asymptotic to the slowest of the two 
defining eigenvectors. 

8) Phase trajectories of systems with one real and two complex 
conjugate roots having initial conditions in the complex eigenplane will 
remain in the complex eigenplane. 

9) Phase trajectories of moderate and heavily damped complex root 
systems, having initial conditions not in the complex eigenplane, spiral 


into the origin of error Space asymptotic to the system eigenvector. 


7S 





10) Phase trajectories of very lightly damped complex root systems, 
having initial conditions not in the complex eigenplane, will spiral 
toward the origin asymptotic to the complex eigenplane. 

11) Phase trajectories of a system with three real roots whose 
initial conditions are not in an eigenplane cannot penetrate an eigen- 
plane and will approach the origin of error space asymptotic to the slow 
eigenvector. 

12) Any third order phase trajectory in response to any initial 
conditions can be described by the superposition of the two phase trajec- 
tories resulting from initial conditions on an eigenvector and initial 
conditions in either a complex or real root eigenplane. (See Chapter IV, 
PP- 67 > 68 ) 

13) If the phase trajectory of a system with one real and a pair 
of complex conjugate roots intersects the eigencone in response to a 
step input, the time of intersection is independent of the size of step 
and the intersection will occur along a faster eigenvector than the 


Original system eigenvector. 
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Chapter VI 
RECOMMENDATIONS 


From the results of this thesis it can be seen that there are 
certain areas of interest which warrant further investigation. A logical 
extension of this thesis would be an attempt to project the results 
obtained herein for third order systems into the more general case of 
Nth order systems. The existence of hypersurfaces in N-dimensional 
error space corresponding to the three-dimensional eigenplanes has 
been well established in the literature. While not yet proven, it 
seems reasonable to assume that there must also exist in N-dimensions 
a hypersurface which corresponds to the eigencone. Thus it should be 
possible to describe the characteristics of phase trajectories of Nth 
order systems in terms of the various hypersurfaces defined by the 


N-dimensional eigenvectors. 


The results of this thesis could also have applications in the 
design of discontinuous control systems. The literature contains a 
great deal of information on various schemes which utilize hypersurfaces 
as switching surfaces, but none of these hypersurfaces are comparable to 
the eigencone. However, in the N-dimensional case, it would have to be 
proven mathematically that under certain conditions the phase trajectories 
intersect an N-dimensional eigencone and that eigencone switching could 


be effected. 


The effect of non-linearities on the phase trajectories of third 
order and higher systems is another field of investigation to which the 
results of this thesis can also be extended. The regions of three- 
dimensional error space in which a system is linear can be defined by 
various planes which establish the boundaries of the linear region. For 
example, in the case of velocity saturation, the linear region is defined 
by the planes E = +kK, where K 1s a constant. Similarly, the linear 
region in the case of torque saturation is established by the planes 
E = +kK. In N-dimensions it seems reasonable to assume that corresponding 
hypersurfaces exist which also establish the boundaries of these linear 
regions, and their effect on phase trajectories has not been adequately 


described. 
276s 
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Appendix B 


ANALOG COMPUTER SOLUTION OF THIRD ORDER DIFFERENTIAL EQUATION 


feneral third order equation: 
Eat pees co =O 


on = -Y, +7 +74 


pen +hG+ GT 


Scaling factors: 
Xe =i 


Ag = KE = Ke = AS 


Basic equation for an integrator: 








Scaling of amplifier stages: 


l. Amplifier 1 
; — 
E=-f{ G@EtbE+cE)dt 
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ow) aX Ee b xe é Va 3 ve . ee 

F=- ; ( 4 « E =n “HE )dt, -[ @é bE+cE )dtp 
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2. Amplifier 3 
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ee Amplifier 5 
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Pipers 2ljeec. co 


No scaling required since these amplifiers are used only as sign changers. 


The analog computer schematic for the general solution is shown in Fig. B-l 
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Appendix C 


COMPUTER SET-UP AND INITIAL CONDITION DERIVATION FOR 
PHASE TRAJECTORIES IN TRANSLATING EIGENPLANE 


The computer set-up as used on the Donner 3100 is given in Fig. C-l. 


The initial conditions for the phase trajectory shown in Fig. 41 


were for a point in the complex eigenplane and were: 
E *=-0.5,E #®+1.0,E = +30 
O O O 
The complex eigenplane was translated to a point on the eigenvector 
corresponding to: 
Eo = + 4,0, E == 2.0; E = 1.0 


Therefore the initial conditions of the point in the translated eigen- 


plane were: 


ie 
li 
7 
+ 
7 
it 


+. 5 


E'SsE+E = +4.0 


These initial conditions gave the phase trajectory of Fig. 42. 
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Appendix D 


DERIVATION OF THE TRANCENDENTAL EQUATION 
cos Bt +o sin Bt = Bs 


~(G4G)t -(r,+r,)¢ 


C ‘a GE 


= : SS D-1) 
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This expression may be rewritten as 
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Dividing through by C, 
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Similarly, 
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Multiplying (D-3) by aK 
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substituting for r,, r r 
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